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Abstract 

Motivated by the evidence^ | 2 | 3 | 4 | 5 | S | in Pr2- 3; Ce a ;Cu04_j / and Nd2-xCe a: Cu04_ 2/ of a magnetic 
quantum critical point at which Neel order is destroyed, we study the evolution with doping of 
the T = quantum phases of the electron doped cuprates. At low doping, there is a metallic 

-t— > 

^ Neel state with small electron Fermi pockets, and this yields a fully gapped d x 2_ v 2 superconductor 

with co-existing Neel order at low temperatures. We analyze the routes by which the spin-rotation 
symmetry can be restored in these metallic and superconducting states. In the metal, the loss 
of Neel order leads to a topologically ordered 'doublon metal' across a deconfined critical point 
with global 0(4) symmetry. In the superconductor, in addition to the conventional spin density 

'— 1 wave transition, we find a variety of unconventional possibilities, including transitions to a nematic 

superconductor and to valence bond supersolids. Measurements of the spin correlation length and 

■^j- of the anomalous dimension of the Neel order by neutron scattering or NMR should discriminate 
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these unconventional transitions from spin density wave theory. 
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FIG. 1: Brillouin zone map showing where the low energy fermions reside for hole-doped [K v ] and 
electron-doped [Q v ] cuprates, as deduced from photo-emission data at low doping. 

I. INTRODUCTION 

Superconductivity in the cuprates emerges on doping an antiferromagnetic insulator with 
either holes or electrons. The hole-doped cuprates generally have higher superconducting 
critical temperatures, but at the same time display a host of complicated phenomena, e.g. 
incommensurate magnetism and charge order, especially in the La series of compounds. The 
electron-doped cuprates on the other hand, provide an interesting contrast, where the phe- 
nomenology appears to be relatively simple. The superconductivity also has (i-wave pairing^, 
but there is no evidence yet for charge order, and the magnetic correlation remain commensu- 
rate even after long-range magnetic order is destroyed. The sharp contrast between electron 
and hole doping must arise from particle-hole asymmetry in Cu-0 planes. The electron-hole 
asymmetry of the Cu-0 plane is evidenced most clearly by photo-emission experiment s 8 * 9 * 10 *^ 
that show a sharp distinction between the Brillouin zone location of the low-energy fermions 
in the very lightly hole- [K v = (±7r/2, ±7r/2)] and electron- [Q v = (tt, 0), (0, n)] doped 
cuprates, see Fig. [Tj 

A further motivation for the study of the electron-doped cuprates is provided by recent 
quantum oscillation evidence for the presence of electron pockets in the hole-doped cuprates 
in a strong magnetic fields It seems most natural to us that these electron pockets reside 
near the Q v . So it seems appropriate to study the physics of the electron pockets where 
they are already present in zero field: in the electron-doped cuprates. Conversely, as we 
will see in this paper, the hole pockets near the K v also play a role in the physics of the 
electron-doped cuprates. Indeed, in both the electron- and hole-doped cuprates, a central 
problem is understanding how the K v hole pockets and the Q v electron pockets reconnect 
to form a large Fermi surface state after the loss of magnetic order. 
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FIG. 2: Schematic phase diagram for the electron doped phase cuprates (following Refs. ll|2|3|4|5|6p . 
The dashed lines indicate finite- T phase transitions. The quantum critical point where Neel (AFM) 
order is lost in the superconductor (SC) is marked with a solid circle. The "doublon metal" is a 
phase proposed in the present paper, which appears when Neel order is lost in the AFM metal; a 
AFM metal/Doublon metal quantum critical point does not appear in the phase diagram above, 
but would be revealed when superconductivity is suppressed e.g. by an applied magnetic field. 
The finite T crossovers can exhibit features of both the AFM+SC/SC and AFM Metal/Doublon 
metal quantum critical points. 



A recent neutron scattering study of the Neel correlation lengtfP in 
Nd2- x Ce x Cu04_j / provides evidence for a quantum critical point at x ~ 0.13, after 
which the Neel correlation length is finite. Remarkably, even at the optimal doping 
x ~ 0.15 (at which long range Neel order is lost) a large Neel correlation length is measured; 
additionally, there is no evidence for incommensurate magnetic order over the entire doping 
range. The relative stability of the commensurate magnetism in the electron doped cuprates 
should be contrasted to the La series of the hole-doped cuprates. In the latter, long range 
magnetic order transforms from the (jr, n) Neel vector to incommensurate ordering vectors 
before being destroyed at dopings typically three times smaller than in the electron-doped 
cuprates. 

These photoemission and neutron scattering measurements suggest the schematic phase 
diagram shown in Fig. [2j as a function of temperature (T) and electron doping (x). 

The focus of this paper is on the nature of the dynamic spin correlations in the electron- 
doped cuprates as a function of increasing doping. It is useful to frame our discussion 
by first recalling the predictions of a conventional spin-density-wave (SDW) theory of the 
evolution of the Fermi surface as function of electron density and the spontaneous Neel 
moment)- 13 * 14 * 15 * 16 ^ . We sketch the results of a mean-field computation in Fig. |3j 



3 



< Increasing AFM order 



i i 



r -i 



i— 1 

1 \ 







n a 






/ \ 







AFM Metal 
or AFM+SC 




Large Fermi 
surface metal 
or SC 



FIG. 3: (color online) Fermi-surface reconstruction at SDW transition showing the presence of 
an intermediate state between the large Fermi surface (rightmost panels) and small Fermi-pocket 
(leftmost panels) states. Following Ref. 15, we have used a band structure appropriate to the 
cuprates, t\ = l,t 2 = 0.32ti and t% = 0.5t2- The right most plots shows the Fermi surface before 
the introduction of a mean-field SDW order parameter. The second from right show the Fermi- 
surface after folding with Asdw = 0, followed by Asdw = 0.05, 0At\ moving left. The top-row has 
chemical potential fi = 0.94, and the bottom row fi = 0.34. The dashed lines indicate the points 
where the d-wave pairing amplitude changes sign in the superconducting state. The AFM+SC 
states in the leftmost panels have fully gapped quasiparticles because the Fermi surfaces do not 
intersect the dashed lines. Similarly, the large Fermi surface SC in the rightmost panels has gap less 
quasiparticle excitations at 4 nodal points, while the intermediate states have 8 nodal points. 



At very low electron doping (x), we have the electron Fermi-pocket states shown in the 
leftmost panels (AFM Metal), with well established Neel order. When this state goes super- 
conducting at low temperature (T), the Fermi surface does not intersect the diagonals along 
which the d x 2_ y 2 pairing amplitude vanishes, and so the resulting o?-wave superconductor 
(AFM+SC) is fully gapped. At large electron doping, we have the large Fermi surfaces 
shown in the rightmost panels, with no Neel order. Now the Fermi surfaces do intersect 
the diagonals at 4 points, and so the d-wave superconductor has 4 nodal points. Examining 
the evolution of the Fermi surfaces between these two limiting cases in Fig. |3j we observe 
that there is generically an intermediate Fermi surface configuration, with Neel order, in 
which the Fermi surfaces intersect the diagonals at the 8 points ±(7r/2, 7r/2) ± (e, e) and 
±(7r/2, — 7r/2) ± (e, — e), for some small non-zero e. The appearance of superconductivity 
at low T will then lead to a <i-wave superconductor with 8 nodal points in the full Bril- 
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FIG. 4: (color online) Analogy between the phases and phase transitions of the physical t-J model 
of electrons and the toy t-J model of S = 1/2 bosons. This figure lists only the metallic phases 



of the electrons, considered in Section III A All non-zero, gauge-invariant condensates (bilinear in 



the g± and z a ) of each phase are noted; those not shown are zero in that phase. The boson-analog 
of the fermionic metallic states are obtained if we replace the g± Fermi pockets by condensates of 



the g± bosons: then the "fermionic Higgs mechanism" discussed in Section III A finds its analog 
in the ordinary Higgs condensate of the g± bosons. Monopoles are suppressed in all phases above. 



louin zone of the square lattice. Thus, in both the metallic and superconducting cases, this 
intermediate state has 8 zero-energy crossings of the fermion dispersion relation along the 
diagonals of the full square lattice Brillouin zone. 

A further motivation for our study is that the 8 diagonal Fermi points of the intermediate 
state are not clearly seen in photoemission experiment d 8 * 9 * 10 ^. Fermi surface crossings are 
seen on only a single point adjacent to the 4 (±7r/2, ±7r/2) points. We therefore explore 
here unconventional routes by which the Neel order at low doping, in the AFM metal and 
the fully gapped <i-wave superconductor, can be destroyed by increasing hole concentration. 

Important aspects of our results on the metallic and superconducting quantum phases 
and phase transitions are summarized in Figs [4] and |5j The right panel of Fig. [4] indicates 
results on a "toy" t-J model of S — 1/2 bosons which we will describe in Section |IV} We 
will see that there is a close analogy between our results for the electronic t-J model and 
the toy boson model, with the latter model having the advantage that duality computations 
of the crossover into confinement can be carried out in explicit detail. 

For the metallic case, we find that the quantum transition out of the Neel state with Fermi- 
pockets (the AFM Metal) is into an exotic 'doublon metal' state without magnetic order (see 
Fig. [4j ; the nomenclature refers to the sites with double occupancy when the Mott insulator 
is doped with electrons. The 'doublon metal' is the particle-hole conjugate of the 'holon 
metal' state described in recent work^^, and both are examples of 'algebraic charge liquids'. 
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FIG. 5: (color online) As in Fig. [4] but for the superconducting phases of the electrons discussed 
in Section |III B| The explicit computations in Section IV are for the boson model, but for paired 
superfluids (which includes all phases above) the results are expected to also apply to electrons. 
The pattern of translation symmetry breaking in the valence bond supersolid is sketched (see also 
Fig. [7]), while that in the nematic superconductor follows Fig. [8^,. Additional forms of translational 
symmetry breaking are also possible in the AFM and non-magnetic phases, and these are discussed 
in Section HV Bl 

These states have topological order and no sharp electron-like quasiparticles. However, they 
are separated from conventional Fermi liquid states by sharp transitions only at T = 0; at 
T > there are only crossovers into the Fermi liquid-like regime. As superconductivity 
always appears as T — > (see Fig. [2]), it is these T > crossovers of the metallic regime 
which are needed for experimental comparisons. We shall show that the spin excitations near 
the transition into the doublon metal are described by a quantum field theory with global 



0(4) symmetry, as indicated in Fig. HI Further, as we discuss below in Section I A, spin 



fluctuations of this 0(4) theory have clear experimental signatures. Section IV will show 



that these metallic phases of the electronic t-J model also have strikingly similar analogs in 
the t-J model of bosons, along with a magnetic ordering transition in the 0(4) class. 

For the superconducting case, we find a number of distinct possibilities, which are illus- 
trated in Fig. |5j From the AFM superconductor we find 3 distinct classes of transitions: 
(a) A transition to a <i-wave superconductor with full square lattice symmetry, which is in 
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the 0(3) universality class. This transition is in the same universality class as conventional 
SDW theory. It is quite remarkable, and novel, that this SDW transition reappears in our 
formalism based on fractionalized degrees of freedom. 

(b) A transition to a d-wave superconductor with co-existing valence bond solid (VBS) order, 
i.e. to a supersolid. Such a valence bond supersolid was initially discussed in Refs. [20][2Tl 
The pattern of the VBS order is columnar or plaquette (see Figs. [5] and [7]), the same as that 
in the insulatoi* 22 * 23 ! (for rational x/2 = p/q with q/2 odd, other patterns of order are possible, 
as discussed in Section IV B 3). This transition is expected to be of the 'deconfmed' variety, 
in the CP 1 universality class, similar to the transition in insulating antiferromagnets^ESl. 

(c) The third transition is described by the CP 1 theory, but with an additional 'doubled 



monopole' perturbation allowed, which will be explained in more detail in Sections IV B 



The non-magnetic superconductor does break the square lattice space group symmetry, and 
the two allowed patterns of symmetry breaking are in Fig. [8j Notice that one them only 
breaks the Z4 rotation symmetry of the square lattice to a Z2 rectangular symmetry, leading 
to a nematic superconductor shown in Fig. [5] 

It is interesting to note that the above possibilities match the patterns of transitions 
found 23 for insulating antiferromagnets as function of the spin S. In particular, case (a) 
occurs for even integer S, case (b) for half-odd-integer S, and case (c) for odd integer S. 
Here we are considering a S — 1/2 antiferromagnet, but with a background of a compressible 
superconductor. As we shall see, the background density fluctuations in the superconductor 
are able to modify the spin Berry phases so that the transitions match those for different S 
in insulators. 

In addition to cases (a), (b), (c), we note briefly that it is also possible that the AFM+SC 
state already has density modulations. Then, the transition involving loss of AFM order 
will lead to modifications in the ordering pattern, as will be discussed in some detail in 
Section IV B An important point is that, in all these cases, the set of allowed periods 
for the density modulations in the supersolid without AFM order are the same as those 
characteristic^ of paired supersolids of density 1 + x. 

An interesting issue, which we shall largely leave open in the present paper, is the nature 
of the spectrum of the fermionic Bogoliubov quasiparticle excitations in the supersolid or 
the nematic superconductor. One natural possibility * 26 * 27 * 28 * is that these excitations initially 
remain fully gapped, as in the AFM+SC state. On the other hand, knowing that the 
supersolid or nematic superconductor has effective density 1 + x, the structure of the Fermi 
surfaces in Fig. [3] suggest that such a d-wave superconductor should have 4 gapless Dirac 
points. In the deconfmed quantum critical theory, the electron spectral function is fully 
gapped along the diagonals of the Brillouin zone. If the gapless nodal points do appear in the 
non-magnetic phase, they would create "Fermi liquid coherence peaks" at the nodal points, 
with the weight of the coherence peak vanishing as we approach the quantum critical point. 
This phenomenon would then resemble that in dynamical mean field theory^EQl, where the 
Fermi liquid coherence peaks of the metal vanish at the metal-insulator transition, revealing 
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a fully gapped single-particle spectrum at the critical point. This issue will be discussed 
further in Section fill Bl 



A. Experimental tests 

We note here that neutron scattering or NMR measurements of the spin excitation spec- 
trum can serve as useful experimental probes of whether the Neel order is lost as in a 
conventional SDW framework, or in a more exotic deconfmed critical point . In particular, 
the temperature dependence of various components of the dynamic structure factor in the 
quantum critical region can measure two crucial exponents characterizing the transition, the 
dynamic critical exponent, z, and the anomalous dimensions of the Neel order parameter, t]n. 
In terms of these exponents, we have^for S N , the zero frequency dynamic structure factor 
at the Neel ordering wavevector (proportional to the elastic neutron scattering cross-section 
at (-7T, 7i")): 

S% ^T { - 2+m)/z ; (1.1) 

for, Sn, the equal-time structure factor at the Neel ordering wavevector (proportional to the 
energy-integrated neutron scattering cross-section at (ir, 7r)): 

S n ~ T { ~ 2+z+riN)/z ; (1.2) 



for £, the Neel correlation length: 

e-T- 1 ^; (1.3) 

The present neutron scattering experimentP only reports the quantum critical behavior of 
the spin correlation length, which is consistent with z — 1. Although data on Sn exists, a 
scaling analysis to extract the exponent Eq. (1.2) has not been carried out. An important 
test of quantum critical scaling would be to check that the exponent that arises from this 
analysis should agree with an extraction of the same index by an analysis of the Cu NMR 
relaxation rate, 

— T m/z . (1.4) 

The values of the exponents in the conventional SDW theory depend upon whether the 
quantum critical region is controlled by a metallic or a superconducting fixed point. For the 
metallic fixed point, we have the Hertz-Millis-Moriya theory^ z = 2 and T]n = 0, while for 
the superconducting case we have the usual 3D 0(3) transition, z = 1 and T]n ~ 0.038. 

Our main new experimentally relevant results in this paper are the values of these ex- 
ponents for the 'deconfined' transition at which Neel order is lost. The exponents depend 
upon whether we are using a superconducting or metallic fixed point, and our results are 
summarized in Table [1} There are no existing numerical results for the 'doubled monopole' 
transition, and so these are not shown: it may well be that this case has a first-order 
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SDW-metal 


SDW-SC 


0(4) 


CP 1 


z 


2 


1 


1 


1 


Vn 





0.038 


1.37 


0.35 



TABLE I: Predictions for the exponents r]N and z by different theories for the quantum critical 
point observed in the electron-doped cuprates such as Nd2_ a; Ce a ;Cu04_ y . Both exponents can 
be measured in experiment by a straightforward analysis of the temperature dependence of the 
equal-time structure factor, as described in the text. The numerical estimates for the anomalous 
dimensions are based on results from previous studies of the 3-dimensional 0(3) [Ref. 33 , 3- 
dimensional 0(4) [Ref. 34 35 and the 3-dimensional CP 1 model inferred from quantum simulations 
of the Neel-VBS transition [Ref. [36] • 

transition. Note the large values of i]n for the deconfined cases, making them clearly dis- 
tinguishable from the SDW cases. In particular, with t]n > 1 for the metallic case, the 
equal-time structure factor, Sn has a singular contribution which decreases with decreasing 
T. 

We also note that for the superconducting case, the properties of the CP 1 field theory are 
not fully settled in the literatur e! 36 * 37 * 38 ^ with a debate on whether the quantum transition 
is second- or first-order. Nevertheless, there is significant evidence 1 22 1 39 1 of a crossover into 
a regime which is described by the CP 1 field theory. Furthermore, even if the transition 
is first-order, it appears to be only very weakly so, and the simulations of Ref. [36] show a 
substantial T > critical scaling regime. 

Because the electron-doped cuprates are always superconducting in the proximity of the 
quantum critical point at low T, the superconducting critical theory described above is the 
correct description at very low-T scales. The normal state theory does however apply at 
temperature scales above the superconducting temperature and hence could be the relevant 
one for experiments over a large temperature scale. An interesting prediction that arises 
from this crossover is that the equal-time structure factor, Sn, could have a non-monotonic 
T dependence. It should first decrease with cooling (when the system is controlled by the 
metallic fixed point with 7]n > 1), and then crossover to increasing with further cooling, 
when the system is controlled by the superconducting fixed point with t]n < 1- 

The outline of the remainder of this paper is as follows: In Sec. [IT] we derive an effective 
field theory for the electron-doped cuprates in a language well suited to discuss both the 
magnetic phases and the non-magnetic ones that appear on the destruction of Neel order. 
In Sec. |III| we discuss the various possibilities for transitions involving loss of Neel order. 
The t-J model of bosons will be introduced in Section [IV] along with a complete duality 
analysis of its phase diagram and its crossover to confining phases. Finally in Sec. |Vj we 
conclude with a summary of our results. 
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II. FIELD THEORY AT LOW DOPING 



We begin with a symmetry-based derivation of a long wavelength effective action for the 
electron-doped cuprates. We will use the low energy excitations of the low doping state to 
build a theory which is valid also at larger doping when spin rotation invariance is restored. 

The motion of a small number of charge carriers in a quantum anti-ferromagnet is usually 
described by the t — J model, 

Ht-j = -J2 M'"';, + h.c.) + J2 J iA -Sj + ... , (2.1) 

where c ia destroys an electron with spin a on the sites i of a square lattice and Si = 
2 c 7 ^a c i/3i with a the Pauli matrices. We shall study the case in which the electrons hop 
on a square lattice. Once extra electrons are doped into the half-filled magnet a constraint 
must also be included. The constraint, 

E c " tc - ^ 1 ( 2 - 2 ) 

is enforced on each site, modeling the large local repulsion between the electrons. It is 
important to note that our results are more general than a particular t-J model, and fol- 



low almost completely from symmetry considerations. The ellipses in Eq. (2.1) indicate 
additional short-range couplings which preserve square lattice symmetry and spin rotation 
invariance. 

Following Ref. HQl but now for the case of electron-doping, we re-write the electron 
operators in a t — J type lattice model in terms of spinons and 'doublons' (for doubly 
occupied sites). Note that here the site occupation is constrained to be J2 a c a c °< — 1- We 
use the following representation for the electron operators, 

c a = e a plF ! g+ (on A) 

c a = -P^ (on B) (2.3) 

where the constraint is b a % a + g + g+ = 1. [We first used c a = e a pb^g on both sub- lattices, 
then rotated the Schwinger bosons on the B sublattice b a — > Sapb , like in the Auerbach- 
Arovas analysis^. Note: e a/3 e^ = —5® 

Now we are in a position to write down the transformation of the lattice fields that 
we have written down under the various square-lattice symmtries and time reversal. We 
require that the composite fields c a transform into each other in the usual way under the 
square lattice symmetries. The implementation of time reversal symmetry is detailed in 
Appendix |X} We thus arrive at the Table [TTJ 

We now proceed to take the continuum limit of the lattice model that we have defined. 
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LATTICE FIELDS: 





T 

- 1 X 


Ddual 
K/2 


rdual 

X 


T 


b a 


e a po 


tP 


E a pO 




T 


e a % 


E a % 


e a % 




9+ 


9- 


9- 


9- 


~9l 


9- 


-9+ 


-9+ 


-9+ 


-gl 



TABLE II: Transformations of the lattice fields under square lattice symmetry operations and time 
reversal. T x : translation by one lattice spacing along the x direction; R^j^- 90° rotation about a 
dual lattice site on the plaquette center (x — > y,y — > —x); I^ ual : reflection about the dual lattice 
y axis (x — ► —x,y — > y); T: time-reversal, defined as a symmetry of the imaginary time path 
integral. The transformations of the Hermitian conjugates are the conjugates of the above, except 
for time-reversal of fermions. For the latter, g± and g\_ are treated as independent Grassman 
numbers and T : g± — > g± . 



In order to do so 23 , we define fields z a = b a + b a and ir a = b a — b\ and integrate out the 
massive Tr a field. We then arrive at the Lagrangian for the z a 

C z = D+z a *D-z a + s\z a \ 2 + u {\z a \ 2 f + ... (2.4) 

where \i = x, y, r is a spacetime index, = ± iA^, A^ is an emergent U(l) gauge field 



linked to the local constraint in Eq. (2.3), and s and u are couplings which can be tuned to 
explore the phase diagram. The Neel order parameter is simply n = z a *a^zp. 

We also need to take the continuum limit for the charge carrying fermions of this model. 
As discussed in detail in Ref. 001 fermions that live on opposite sub-lattices carry opposite 
charges under the gauge field, A^, and hence must be represented by two distinct continuum 
fields g± (both fields are centered at the lattice momentum, Qx). The lowest derivative term 
consistent with the symmetry of the g± is, 

^E^-M-^V (2-5) 
where m is the curvature of the fermion bands and q = —q. Finally, by requiring consistency 



with the lattice transformation properties of the continuum fields, presented in Table III 
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CONTINUUM FIELDS: 





T 

- 1 X 


E>dual 
H */2 


rdual 

X 


T 


z a 






e a pz p * 


E a pzP* 


9+ 


-9- 


-9- 


-9- 


-A 


9- 


9+ 


-9+ 


9+ 


-9* 



TABLE III: Transformation properties under square lattice symmetries and time reversal, of con- 
tinuum fields entering the effective action. Conjugate fields transform into the conjugate of the 
transformed fields except for T : g± —* g± 

the lowest allowed derivative term that couples the opposite fermions g± can be deduced, 



C z - 9 = A e af3 



+ £a(3 



g\ (D+g_) z a (D-z p ) - g\ (D+g_) z a (D~z p )] (2.6) 
'gl (D-g + ) z a * (D}z?*) - gl (D~ g+) z a * {D+z^*)} + c.c. 



This is the analog for electron-doped cuprates, of the well-known Shraiman-Siggia termP^ 
in the hole-doped case. Remarkably, this term has two spatial derivatives; there is no term 
allowed with a single spatial derivative (as is found from a similar analysis in the hole-doped 
case^ see alscP^). The extra derivative makes the effect of this term weaker. The weakness 
of this coupling, which arises because of the BZ location of the low energy fermions, (which 
in turn is ultimately tied to the p-h asymmetry in the Cu-0 layers) is the fundamental reason 
for the robustness of the commensurate Neel correlations in the electron-doped cuprates as 
compared to the hole-doped case. These correlations extend at least up to optimal doping^ 
and possibly beyond giving us confidence in the present approach. 

The complete effective action for the electron-doped antiferromagnet is then S = 
J d 2 rdr(C z +C g +C z ^ g )+SB- The final term, Sb contains the Berry phases of the monopoles, 
and has the form 

<S B = ^X>^ (2-7) 

3 

for monopoles with integer charges rrij on the sites j on the dual lattice; (j is fixed at 
C = 0, 1, 2, 3 on the four dual sub-latticed^. 



A. Neel order and superconductivity 

We now discuss the phase diagram of the field theory presented in the previous section. 
Some of the analysis parallels that presented in Refs. SU] and HH] for the hole-doped case. 

The phases are most easily characterized by using a representation for the physical elec- 
tron annihilation operator \P a (r) in terms of the fields we have introduced above. We first 
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express the electron operator in terms of its components at momenta at Q\ and Q2, 



* a (r) = e l ^ la {f) + e i ^ 2a (r). (2.8) 

Then, as in Ref. HOI we use the symmetry transformation properties to deduce the unique 
bilinear combination of the fermion and and CP 1 fields that transform in the way that the 
physical electrons ^2 should, 

*i,2a = e a pz p *g + Tz a g-. (2.9) 

The phases is Fig. [2] can now be characterized in terms of the z a and g± degrees of free- 
dom: 

(i) AFM metal : This is the Higgs phase of the gauge theory, in which there is Higgs con- 
densate of z a with (z a ) 7^ 0. As discussed in Ref. HQJ the "Meissner" effect associated with 
this Higgs condensate ties the gauge charge to the spin quantum number. So for Neel 
order oriented along the z axis, the g± fermions carry spin S z = ±1/2 and reside in Fermi 
pockets. The resulting phase is then identical to the AFM metal phase obtained in SDW 
theory, and shown in the left panel of Fig. |3j 

(m) Doublon metal : This is the particle-hole conjugate of the holon metal, and is a non- 
Fermi liquid 'algebraic charge liquid'. We have (z a ) = 0, and and the phase is described by 
the gapped z a quanta and the g± Fermi pockets interacting via exchange of the gauge 



force. We observe from Eq. (2.9) that the physical electron involves a convolution of the 
propagators of the z a and g± , and so will not have Fermi liquid form. We note that the holons 
and the holon-spinon bound states, discussed in previous work^^l on the holon metal, are 
also legitimate excitations of the doublon metal. Here they are likely to be gapped, but at 
T > will contribute photoemission spectral weig htP near the K v points in Fig [TJ 
(in) SC phases: As discussed in Ref. [TSJ the nearest-neighbor hopping term, and the gauge 
forces, will induce a pairing of the g± fermions. Let us assume a pairing of the form 

(g +1 (k)g_ 1 (-k))=A(k). (2.10) 



Then the pairing signature of the electrons can be computed from Eq. (2.9) and (2.10): 
the various possibilities are discussed below. If we also have (z a ) 7^ 0, then we obtain the 
AFM+SC phase of Fig. [2} This is a stable phase, because the z a Higgs condensate quenches 
the gauge fluctuations and also the monopoles; its physical properties are identical to the 
AFM+SC phase obtained in the SDW theory noted in Fig. [3} A superconducting phase with 
(z a ) = is the doublon superconductor, and this is not stable: proliferation of monopoles 



will lead to confinement, as we shall discuss in Section [HI B 

The remainder of this subsection will characterize the symmetry properties of the possible 
superconducting phases. We also allow long-range Neel order by a condensation of the 
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d xy d x 2_ y2 g 




FIG. 6: (color online) Various order parameters that are even under inversion A(fc) = A(— k). 
Shaded areas are positive and white areas negative; Thick red lines denote zeros of the order 
parameter. The first (second) row has order parameters that are even (odd) under translation by 
a Neel vector A(k + if at) = — A(fc). The columns are order paramteres that transform similarly 
under square lattice operations: A(/c) — > ±A(/c) under rotation by ir/2 and reflection across x,y 
axis. Only the second row can appear in the co-existense phase discussed in the text, in which 
pairing is between electrons on opposite sub-lattices. 



CP 1 fields with (z a *zp) = 5 a pm a : the Neel order is then polarized in the z direction with 



spontaneous moment m-f — m^. Using these averages, Eq. (2.10), and the expressions for the 



physical electron operators, Eq. (2.9), we can calculate the required anomalous averages, 





-AO) 




-*)> 




-*)> 



£ a p[mpk{k) + m a A(-k)\ (2.11) 

At the critical point from the AFM+SC state to the SC state the Neel order parameter 
vanishes, i.e., = and the correlator should disappear (this follows from the 

restoration of full translational invariance), indicating that A(k) = A(—k). Since the super- 
conducting instability arises out of a short-range attractive interaction it is most natural to 
expect s-wave pairing. Remarkably, this naturally leads to d x 2_ y 2 pairing for the physical 



electrons [as can be verified from Eq. (2.11) by substituting A(k) = Aq]. However since the 
underlying g± particles are in an s-wave state, the quasi-particles in this d x 2_ y 2 supercon- 
ductor are fully gapped. We propose that this is the quantum state that describes the phase 
AFM+SC in Fig. [2]and that is observed in the region of co-existence in Nd 2 _ I .Ce a; Cu04_ 3/ ^. 
We note that with increasing x, the Neel order is suppressed making the gauge field medi- 
ated attraction (that causes superconductivity) stronger, which in turn is expected to result 
in an increase of T c , consistent with experimental observations. For the sake of complete- 
ness, we present the other symmetry allowed options for pairings (see Fig. [6] second row): 
A(k) = k 2 x — ky corresponds to the s case, A(k) = k x k y {k 2 x — k y ) corresponds to the d xy case 
and A(k) = k x k y corresponds to the g case; all these states have nodal excitations. Finally, 
the condition A(k + K]y) = —A(k) that is satisfied by all the order parameters deduced from 
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Eq. 2.11 (and that are illustrated in the second row of Fig. [6]) follows quite simply from the 
fact that the phase factor e lKN ' r is +1(— 1) on the A(B) sublattices and that pairing occurs 
only between electrons on opposite sub-lattices. 



III. QUANTUM CRITIC ALITY 

We now turn to our main results on the quantum phase transitions involving loss of 
Neel order as described by the low energy theory introduced in Section [TTJ the results were 
summarized in Fig. [4] and [5] 



A. Metallic states 

First, let us consider the transition without superconductivity, destroying magnetic order 
in the AFM metal, leading to the doublon metal. This transition is described directly by 
the field theory in Section [XXJ and is associated with the condensation of the z a spinons in 
the presence of the g± Fermi surfaces. 

At T = 0, such a transition between metallic states could be induced by destroying super- 
conductivity by an applied magnetic field. Moreover, even at zero magnetic field, the quan- 
tum critical region at temperatures above the superconducting T c could be controlled by the 
crossovers of an underlying AFM metal/doublon metal quantum critical point. Monopoles 
can be ignored in the following analysis because they are suppressed by the gapless excita- 
tions at the g± Fermi surfaces^ (see Appendix [b]) . The resulting state without antiferro- 
magnetism therefore carries gapless gauge excitations, and as we noted earlier, realizes an 
algebraic charge liquid which we call a doublon metal. 

The theory for this transition follows the analysis of a formally similar transition of bosons 
and fermions coupled to a U(l) gauge field in Ref. HE, In this previous case, the bosons were 
spinless and fermions carried spin, whereas here the fermions are spinless while the bosons 
carry spin. However, for the quantum criticality, the more significant difference is that the 
quadratic action for the z a bosons has a relativistic structure, unlike the z = 2 dispersion 
in Ref. 051 

The renormalized gauge field propagator is a key ingredient in our analysis. This 
depends upon the polarizabilities of the g fermions and the z a bosons at the quantum 
critical point. We evaluate these from the bare actions C g and C z , and will confirm later 
that the same results hold in the fully renormalized critical theory. As usual, the fermion 
polarizability screens the longitudinal A^ fluctuations, and the only potential singularity 
arises from the transverse A^ propagator, D. In the Coulomb gauge, this has the low 
momentum and imaginary frequency forrrP^ 
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Here the \oo\/k term in the denominator is the contribution of the g fermions, and is present 
for \u\ < vpk, where vf is the Fermi velocity. The k term emerges from the critical z a 
correlator (it coefficient is proportional to the critical conductivity of the z a 's). 



Let us now compute the consequence of the overdamped gauge fluctuations in Eq. (3.1) 
on the z spectral function. At leading order the z a self energy at criticality is 

Z z (p,ie) ~ [du> I £ k ^-^-K?lp 1 

J J k + X \uj\/k (uj + e) 2 + (k + p) 2 v ' 

It is now easy to confirm that this expression for Z z is non-singular at low p and e, and does 
not modify the leading behavior of the z propagator. In particular, the on-shell self energy 
has the imaginary part 

lmS z (p,e = p) ~ p 3 , (3.3) 

which is clearly unimportant to the critical theory. Thus the overdamping of the gauge 
fluctuations by the g± fermions strongly suppresses their influence on the z a excitations. 



Indeed, in the z = 1 scaling, the k term in the denominator of Eq. (3.1) can be neglected, 
and the renormalized action for the transverse component of the gauge field is ~ Ay(|ct?|/fc); 
this scales as an anisotropic "mass" term for the gauge boson. Thus we can view this feature 
as a fermionic version of the Higgs mechanism, in which the low energy excitations of a Fermi 
surface quench the gauge field fluctuations. We will comment further on this analogy with 



the Higgs mechanism in Section IV 



In a recent work 47 in a different context, Senthil has computed the consequences of the 



singular interactions associated with Eqs. (3.1) and (3.2) on the spectral function of the 



fermions, and the associated formation of critical Fermi surfaces^. All those results apply 
here too to our theory of the transition from the AF metal to the doublon metal. 

At this point, we are now prepared to integrate out the gauge boson and the g± 
fermions, and obtain an effective theory for the z a spinons. Keeping only the terms poten- 
tially relevant near the critical point, the resulting effective action has the structure 



* ff = / d 2 rdr 



\d^z a \ 2 + s\z a \ 2 + u{\z a \ 2 ) 2 \ + A J d 2 kdu[\z a \ 2 ](-k,-Lu) M [\zp\ 2 ](k,co) 

(3.4) 

The last A term is a consequence of the compressible fluctuations of the g± Fermi surfaces, 
which couple to \z a \ 2 via a contact ternP. At A = it is now evident that S z s describes 
a transition for the loss of Neel order by a conformal field theory in the 0(4) universality 
class. We can therefore ask for the scaling dimension of A at this conformal critical point. 
This follows from a simple scaling argument^ 

dim[A] = -3 + - (3.5) 



The 0(4) model has^" v = 0.733 and so A is an irrelevant perturbation. Further, when 
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we account for the long-range Coulomb interactions between the g± fermions, there is an 
additional factor of k in the A term, and A is then more strongly irrelevant. 

We have now established that the transition from the Neel-ordered Fermi-pocket metal 
to the doublon metal is in the 0(4) universality class. The Neel order parameter itself is a 
quadratic composite of the z a . It transforms under the symmetric, traceless, second-rank 
tensor representation of 0(4), and the scaling dimension of this composite operator has 
been computed earlier^^l. From the field-theoretic analysis of Calabrese et al. 34 we find 
fjfq = 5 — 2y 2 .2 — 1.374(12), while the Monte Carlo simulations of Isakov et al¥^ we obtain 
r] N = 1.373(3). 



B. Superconducting states 

Now we discuss the transition out of the superconducting AFM+SC state with increasing 
doping. Because the g± fermions are fully gapped in the superconductor, they can initially 
be ignored in the analysis of the critical theory. The remaining z a excitations are described 
by the CP 1 model. So a natural initial guess is that the critical theory for the loss of Neel 
order is the same as that in the insulator^!. The presence of superconductivity here does 
induce additional gap less density fluctuations, but these are irrelevant!^ as long as v > 2/3 
for neutral systems, and generically unimportant with long-range Coulomb interactions. 
Further, the paramagnetic state so obtained is not a BCS superconductor, because the 
fermionic Bogoliubov quasiparticles carry no spin. Rather, as discussed in Ref. (TSJ it is 
a "doublon superconductor". However, once we have moved away from the critical point, 
there are no gapless excitations which can serve to suppress monopoles in the U(l) gauge 
field. We expect that the condensation of the monopoles at a large secondary length scale 
will induce confinement, leading to a generic instability of the doublon superconductor. 

We are interested in the nature of the confined state. For the corresponding transition in 
the insulator^!, the confining state was the valence bond solid (VBS) which was induced by 
the Berry phases on the monopoles. However, here there is the possibility that the density 
fluctuations of the superconductor can modify the influence of the Berry phases. Because 
the g± fermions are paired in the AFM+SC state and at the critical point, it is plausible that 
a t-J model in which the g± are bosons (and the z a remain bosons) should have essentially 
the same properties in its charge and density correlations in their respective paired states: 
we are merely replacing the internal constituents of the Cooper pairs, but this should not 
modify the nature of the phase and vortex fluctuations of the superfluid. We will examine 



such a t-J model of bosons in Section IV we are able to carry out an analysis of the influence 
of monopole condensation in some detail, and find 3 distinct possibilities which were listed 
earlier in Section |T] and Fig. |5] 

(a) A conventional 0(3) transition, as in SDW, theory, to a <i-wave superconductor with 
full square lattice symmetry. The monopole Berry phases are precisely cancelled by density 
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fluctuations in the superfluid, and so the monopoles confine the z a spinons into the vector 
SDW order parameter. 

(b) A 'deconfined' CP 1 transition to a valence bond supersoli d 20 * 21 ^, where the pattern of the 
VBS order is the same as that in the insulato r^ 2 * 23 ! (see Fig. [7]). Here the monopole Berry 



phases remain as in Eq. (2.7). For rational x/2 = p/q with q/2 odd, other patterns of order 



are possible, as discussed in Section IVB3| 



(c) A direct transition to a <i-wave superconductor with square lattice symmetry broken as 
in the states in Fig. [8j one of which is a nematic superconductor. In this case, the monopole 
Berry phases are only partially compensated by the superfluid modes, so that monopoles 
with even magnetic charge are allowed at the transition. Little is known about the critical 
properties of such a 'doubled monopole' theory, and it is possible the transition is first order. 

The above list exhausts the possible transition out of the AFM+SC state, for the case in 
which the AFM+SC state does not have density modulations of its own. In Section IV, we 
consider the further possibility that the AFM+SC state already carries density modulations 
(so that it is also a supersolid). We classify transition out of such states: monopole conden- 
sation modifies the nature of the density modulations in the non-magnetic supersolid, and 



these will be described in Section |IVB| An interesting feature of the resulting supersolids 
is that the set of allowed wavevectors for density modulations are the same as those of a 
model of paired particle j 25 * 51 ! of density 1+x. Thus, once Neel order is lost, the primary role 
of the monopoles is to account for the 'background' density of one particle per site in the 
Mott insulator, and to combine this density with the doped particles to yield states which 
are sensitive to the total density. 

This sensitivity to the total density in the supersolid bears some similarity to the con- 
straints placed by Luttinger's theorem on the volume enclosed by the Fermi surface in Fermi 
liquid states®'. Let us explain the connection more explicitly. We can view the monopole 
Berry phases as arising from a filled band of anti-holons in the insulator, and these are 
extracted into the confining electronic states^. To see this, let us recall the origin of the 
monopole Berry phase in Eq. (2.7). This can be traced to the constraint b]j) a -\- g + g+ = 1 ap- 



plied on every site of sublattice A (there are parallel considerations on sublattice B, which we 
will not write down explicitly), and implemented by a Lagrange multiplier A in the effective 
action with the term 

' dr\(bib a + g{g + -l). (3.6) 



The fluctuations of A are A T on sublattice A (— A T on sublattice B), and the —1 in the 



brackets above evaluates^ to Eq. (2.7) for a monopole configuration of A^. Let us now also 



allow for the gapped holon states^^ by the holon operators f±, in which case the term in 



Eq. (3.6) generalizes to 

i J dr\(bib a + gig + + flf + -l). (3.7) 
Finally, we perform a particle-hole transformation to anti-holons h + = /| (which are distinct 
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from the doublons) to obtain 

i J dr \(b]j> a + g\g + - h\h + ). (3.8) 

In this form, there is no —1 in the bracket, but we have a filled band of h + anti-holons — thus 
we have an alternative book-keeping in which there is no monopole Berry phase, but we do 
have to account for the unit density of anti-holons in the Mott insulator. After confinement 
with spinons, it is this density which contributes to the expansion to the large Fermi surface 
in the Fermi liquid, and the sensitivity of the supersolid to density 1 +x. 

Next, we consider the structure of the physical electron spectral function in the supersolid. 
We focus on momenta along the diagonals of the square lattice Brillouin zone. Right at the 
critical point, the g± and f± fermionic excitations are fully gapped, and so the electron 
spectral function (which is a convolution of these fermion Green's functions with those of 
the z a ) is also fully gapped. Moving on the confining side of the critical point, a natural 
possibility is that electron spectrum remains fully gapped 2 - 6 * 27 * 28 * . However, given the fact 
that the confining supersolid consists of a density 1 + x (as in a 'large' Fermi surface), from 
Fig. [3] it would seem natural that this state has 4 gapless nodal quasiparticles; if so, the 
total spectral weight in these low energy fermions would vanish as we approach the critical 
point, in a manner we expect is related to the scaling dimension of the monopole operators. 
It is interesting to note that this vanishing of low energy fermionic spectral weight resembles 
the phenomenon of spectral weight transfer in dynamical mean field theory 2 ™. We also 
note that the emergence of gapless composite fermions from gapful constituents (here the 
'composite' electrons consist of gapped spinons and holons) has counterparts in a number 
of particle theory modekP^l 

It is useful to discuss this theory in the context of recent ideas by SenthiP^ on 'critical 
Fermi surfaces'. In the latter framework for a transition to a o?-wave superconductor with 4 
nodal points, the nodal fermions would be part of the critical theory, and then the deconfined 
critical theory would not be the CP 1 model. Such a scenario would be realized here if the f± 
holon Fermi surfaces formed in the AFM state (this is compatible with current photoemission 
experiments^), and the magnetic disordering transition led to a holon superconductor with 
gapless Dirac excitations^. A confinement transition on the holon superconductor would 
then realize this scenario, but only at 'critical Fermi points' and not on a 'critical Fermi 
surface'. 



IV. t-J MODEL OF BOSONS 



The Higgs-like suppression of the fluctuations in Section III A suggested to us that we 



examine a toy model of bosons obeying the same t-J model described here. In other words, 
we will consider the same theory presented in Section [TT] but now the g± and the \l/ fields are 



19 



all bosons. We can make quite reliable statements about the phases of this model, including 
the role of monopoles and Berry phases. 



A further motivation to examining this model was noted in Section III B this is an 
efficient way to analyze the paired superconducting states, where we expect pairs of bosons 
or fermions to have similar properties. 

The analogy between the phases of the electronic model and the toy boson model were 
summarized in Fig. [4] and [5] The parallel of the Higgs-like effects in the metallic phases 



of Section III A appears when we replace the g± Fermi pockets by Bose condensates of the 
g± — the corresponding transitions in the boson model are then in the same universality 
class as in the metallic electronic model. As indicated in Fig. |5j the boson model also has 
parallels to the transitions of the superconducting sector of the electronic model which were 
discussed in Sections IIII B[ This will be described in more detail below. 

First, let us list the phases of the boson t-J model of interest to us: 
(i) AFM boson superfluid: Here both the z a and the g± condense with 

(z a )?o , (g±)^o. (4.1) 

The presence of these condensates implies that both and monopole fluctuations are 



suppressed, as in the AFM metal. Also, by Eq. (2.9), the physical boson operator \l/ also 
has a non-zero condensate. So this state has AFM order and a flux quantum of h/e. 
(ii) Paired boson superfluid: Now spin rotation invariance is restored, with 

(z a ) = , (<7±)^0. (4.2) 

However, the g± condensate is sufficient to continue to suppress both the and the 
monopole fluctuations, making this state the analog of the doublon metal. Two other 
characteristics of this state reinforce the analogy with the doublon metal: (?) the z a quanta 
represent stable, neutral, S = 1/2, gapped excitations, which are also found in the dou- 
blon metal, and (ii) the action of an isolated monopole diverges linearly with system size 
because of the Higgs condensate, and a similar linear divergence appears^ in an RPA-like 
estimation^ of the monopole action in the doublon metal (see Appendix [b]) . With the con- 



densates as in Eq. (4.2), as discussed in Ref. EH the only gauge-invariant condensate carries 



charge 2e, and so the flux quantum is h/(2e). A comparison of Eqs. (4.1) and (4.2) shows 
that the transition between the AFM boson superfluid and the paired boson superfluid in- 
volves criticality of z a alone. The A^ mode can be ignored and so it is evident that the 



critical theory is the 0(4) model in Eq. (3.4), but with the last density fluctuation term 



replaced by the analogous term for a superfluid^. The latter term is also irrelevant, by an 
argument similar to that made for the electronic case. 

(Hi) AFM paired boson superfluid: Now we condense the z a , but only allow for a paired 
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condensate of the g± bosons with 



(z a )^0 , (g±)=0 , (g + g-)^0. (4.3) 

There is antiferro magnetic order, and the flux quantum is h/(2e). The z a condensate is 
sufficient to suppress both the and the monopole fluctuations, making this state the 
analog of the AFM superconductor in the electronic model. In some cases (to be discussed 



below in Section IV B) this state will also break translational symmetry, i.e. it will become 
a supersolid. 

(iv) Paired boson supersolid: The only condensate is that associated with the paired bosons: 

(*«) = , (g±)=0 , (g+g-)^0. (4.4) 



This is the most interesting state here: the and monopole fluctuations are not suppressed, 
and we expect a crossover to a confining state. The same phenomenon also appeared in 
the electronic case with the doublon superconductor, which we argued was unstable to 
confinement to a <i-wave superconductor. The key advantage of the toy boson model is 
that we can describe the crossover to confinement in some detail, as will be presented in the 
following subsections. Our main result will be that there are generally periodic bond/density 
modulations in this phase, i.e. it is a supersolid; we include here the case of the nematic 
superconductor, in which only the rotational symmetry of the square lattice is broken. 
Finally, we will demonstrate that these modulations are characteristic^ 1 ^ of the total density 
of bosons, 1 + x. 



A. Duality and symmetry analysis 

We will apply the analog of the duality methods presented in Refs. l2~4B25ii51l to this model. 
These dualities are only operative for abelian symmetry, and so we shall replace the SU(2) 
spin symmetry by a U(l) symmetry of spin rotations about the z axis. 



We write the spinons, z a , and represent them by two angular degrees of freedom = e 



z l = e%di - Similarly we take the g± (which are now bosons) and write them as g± = e t ^ >± . 
These fields are coupled to a compact U(l) gauge field A^, with the same charges as in 
the body of the paper. Finally, the monopoles in A^ are endowed with the Haldane Berry 
phases^^in Eq. (2.7), to properly include the physics of the insulating antiferromagnet. 



The simplest model consistent with such a framework is written below. Here we have 
discretized spacetime onto the sites of direct cubic lattice with sites j and is a discrete 
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lattice derivative. 



Z = Y\ J d9yd9^jd(p + jd(f)_jdAj^ exp | 



+ -J X] COS ( 6 A«f aAj/Aj-a) - 5 B j . 
6 □ / 



(4.5) 



Apart from the coupling constants, fT, L, e 2 , the action contains two fixed external fields. 
The uniform static external electromagnetic field = ijld^r, where Jl is the chemical 
potential; the value of /Z is adjusted so that density of each g± boson species is x/2. The 



last term accounts for the Berry phases linked to the monopoles in An by Eq. (2.7). 



To be complete, we should also add to Eq. (4.5) a staggered chemical potential which 
preferentially locates the g± on opposite sublattices, as has been done in previous work^^. 
However, this term is not essential for our conclusions here, and so we omit it in the interests 
of simplicity. 



The duality analysis of Eq. (4.5) is most transparent when the action is written in a 
Villain (periodic Gaussian) form. We do this by introducing the integer- valued fields, p-^, 

PijiJ-i ^+jM> 

which reside on the links of the direct lattice, and the integer-valued q m 
which resides on the links of the dual lattice. The dual lattice sites are labeled by j. 

z = Y Y Y Y EII /^^ii#+j#-i% ex p| 

{PTjmI {PUi*} WjaJ {"-jm) 3 

J7 Y ( A ^+i - A 3» ~ B m ~ 2v™ +jAt ) 2 ~^jY ( A /^-i + A J> _ 5 J> _ 



- i Y (^ A « A 3* - 2 ^) 2 -jY J • ( 4 - 6 ) 

□ J / 

An advantage of this periodic Gaussian form is that we are able to write an explicit expression 
for the monopole Berry phased the fixed field C, = 0, 1, 2, 3 is the same as that appearing 



in Eq. (2.7) 



Now we proceed with a standard duality transformation of this action. Initially, this maps 
the theory onto the integer valued spin currents Jy^ and Jy^, the integer valued charge 
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currents H + ja and H-ju, and the integer valued fluxes Q m with the partition function 



Z d = J2 E E E E Constraints eX P I ^ + J lm) 

{•W {-W {H-j„} {H- m } {Q m } \ m 

L 2 / 1 \ 2 \ 

- 2 E + - y E ( ^ - 4 A ^ J - * E s * + ( 4 - 7 ) 

j> w in I 

The summations in Z& are restricted to integer- valued fields which obey the local constraints 

A„ Jti„. = , A„Ju-,, = , A, .Ha- ,, = , A„iy_ 1 „. = , 



£[iu\A. u Qj\ — J^jn + Jjj^ + — H_j^. (4.8) 

We solve these constraints by introducing the dual gauge fields and a; J(U whose fluxes are 
the spin currents, the dual gauge fields b +m and 6_ JAt whose fluxes are the charge currents, 
and a height field h 3 whose gradients are the fluxes. Finally, we promote these dual 
discrete fields to continuous fields by introducing the dual matter fields e~ ia ^ and e~ ta ^ 
which annihilate vortices in the z^i spinons, the dual matter fields e~ lf3+J and e~ tf3 ~ J which 
annihilate vortices in the g± charged bosons, and the corresponding vortex and monopole 
fugacities. This leads to the dual theory in its unconstrained form 

Z d2 = Y\ J da^^daij^db+j^db^j^dhjda^da^dp+jdp^j exp | 

- vE ((^aA^a) 2 + (c^aA^a) 2 ) 
A □ 

r / / — \ 2 / _ X 2 N 

L sr-^ If . /i \ / jJL 



~~ \ E f ^aA^+ja - + ( 6 ^aA 

m 
m 

+ Vvc ^ ( C ° S ^^+3 ~ 27lb +3v) + COS _ 2nb -3^)) 

3V 

+ y m J2 cos ( 27r ^ + a Tj + «Ij + P+3 - P-3 + |Cj) j • 
The average flux of 6± is /Z/L and this should equal half the electron density, x/2. 



(4.9) 



The action in Eq. (4.9) appears to be of daunting complexity, but its physical interpre- 



tation is transparently related to the direct theory. There are 4 vortex matter fields, e"* T , 
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TABLE IV: Correspondence between local operators in direct theory and monopole operators in 
dual theory. q a , q^ and q c are monopole fluxes assoicated with gauge fields a, b and c respectively. 
The subscript s labels spin (s = 1 for j and s = — 1 for j) 



e"* 1 , e i/3+ , e l/3_ . These annihilate vortices in z^, g + and g_ respectively. These 4 matter 
fields carry unit charges under 4 U(l) gauge fields, a-f, aj_, &+, and 6_ respectively. Of these 
4 gauge fields, one combination is always Higgsed out by the scalar field h (by 'Higgsed' we 
mean that the gauge boson acquires a mass via the Higgs mechanism) . The latter is related 
to the monopole annihilation operator e 2mh , and the monopoles carry Berry phases e m ^l 2 . 



Let us now make a further simplification of the dual action in Eq. (4.9). As the gauge 
field combination a-f + aj_ + b + — 6_ is always Higgsed by the h field it is convenient to 
integrate these two fields out, obtaining the dual action 

S d = ir^(e^ A A ! ,a A ) 2 + L^(e^ A A,6 A -|(5 / , T ) 2 + (ir + L)^( e ^ A A 1/ c A ) 2 
□ □ □ 

- Vvs (cos(A At a b - 2ixa 3il - 2ixc m ) + cos(A A1 a ij + 2ixa ni - 2irc 3ft )) 

- Vvc ^2 ( cos ( A m^+j ~ 2r{h m + 271 " c jm) + cos(A M /3_j - 2nb m - 2nc m )) 
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- y m cos(a Tj + a l3 + (3+ 3 - (3. 3 - -Q (4.10) 

3 

The resulting action has three gauge fields: a, b and c. The flux of a is related to the 
magnetization density, S z = e TU \A u a\, the flux of b to the electron pair density, n/2 = 
e TU \A u b\. Finally the field c introduces interactions between spinon and doublon vortices. 
When c is Higgsed out (as happens in the paired boson superfluids), the e ll3± are the physical 
vortices in the superconducting order parameter^ which carry flux h/(2e); otherwise they 
carry flux h/e. As usual, gauge invariant local operators in the direct picture correspond to 



monopole operators of the dual gauge fields (see Table IV ) 



For notational convenience below, we define the spinon vortices ip^i by 

^ T = e ia T , ^ l = e ia K (4.11) 
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1. Symmetries 



A crucial part of our analysis will be an understanding of the symmetries of the action 



in Eq. (4.10) 



First let us consider the action of the space group symmetry of the square lattice. Fol- 
lowing the analyses of Refs. |25j EH we will consider the operations T X)V (translation by one 
lattice site in the x,y directions), R^/f (rotation by a 90° about a dual lattice site), J^? ual 
(reflection x — > —x, with the origin on a dual lattice site), and T (time-reversal). The action 
in Eq. ( |4.10 ) is invariant under these operations with the transformations: 



T x : ip-^ — > iip\, ip l -> 
T y : ^ T ^} 



The non-trivial transformations of the spinon vortices above are a consequences of the 



monopole Berry phases, Q in Eq. (4.10). 

We will be interested below in taking the continuum limit of the effective action for these 
fields. Here we have to be careful about the fate of the Cooper pair/doublon vortex fields 
e tf3± . Indeed, our vortex fields e l/3± are propagating in the background of an average flux for 
the b field that is dual to a finite electron density x/2. We will work at a rational density, 

K <" 3 > 



where p and q are relatively prime integers, and then (as discussed at length in Ref. 
there are q degenerate minima in the Hofstadter dispersion. We label the vortex excitations 
at these minima by the complex fields <p±i, with I = 0, 1, 2, ...q — 1. Thus, in the continuum 
limit, the vortex fields e l/3± are replaced by the 2q fields ip±i. Moreover, once the fields 
e »p± split into <p±i multiplets, the transformations become even more non-trivial due to the 
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presence of a background flux of the b field, 



T x : 
T ■ 

ly . 



r>dual 
K n/2 ■ 

rdual . 

x 

T: 



(fal — > UJ~ l (fal 
1 

v<7 m=0 



(4.14) 



where all indices are implicitly determined modulo q, the "bar" operation exchanges + 
and, 



uj = e 



2-nip/q 



Notice that for the transformations in Eq. (4.14) we have 



T T 



UT X Ty 



(4.15) 



(4.16) 



and this algebra is cruciaP^ in ensuring the g-fold degeneracy of the vortex states. The 
factor ui is understood as a transformation in the U(l)b symmetry group (see below). 

In addition to the space group operations, we should also consider the symmetries asso- 
ciated with global parts of the 3 U(l) gauge groups. More explicitly, we define the transfor- 
mations in the global parts of the gauge groups as, 



<p +l -> e* e (p +h cp_i 



V; 



e<Vi 



e v 9 - 



(4.17) 
(4.18) 

(4.19) 



Then, combining the transformations of spinon vortices Eq. (4.12) and Cooper pair/doublon 

(4.20) 



vortices Eq. (4.14) 



T T 

± y ± x 



-l) c (-uj) b T x Ty 



with all the other relations in the lattice group as in the non-projective case. 

We have now enumerated all the symmetries which will determine the structure of the 
effective action and the phases. However, these symmetries are still somewhat cumbersome, 
and it is useful now to define certain bilinears whose transformation properties are somewhat 
simpler. 

First, we define bilinears of the Cooper pair/doublon vortices. We introduce a set of pair 
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vortex operators^, 



9-1 

7m „ = U mn/2 V-lV+J+n 



US' 



I III 



with the transformation properties, 



1=0 



(4.21) 



T 

r>dual 
K n/2 

rdual 

x 

T 



^Ymn 

linn 
Iran 
'"Ymn 



!—n,m 
7— m,n 
7mti 



'fmn 
'fmn 
'"Ymn 



^Ymn 
^Ymn 

e -2ie 



(4.22) 



Note that the space group transformations of the ■y mn are just those of the physical particle 
density operator at the wavevector (27rp/q)(m,n). However, a crucial point is that j mn 
is not equivalent to the particle density operator: this is a consequence of the non-trivial 
transformation of 7 mn under U(l) c above. Only combinations which are neutral under U(l) c 
are physically observable. 

Next, we consider the following bilinear of the spinon vortices 



(4.23) 



which has the transformations 



T 

± X 

T 

rjdual 
rdual 



T 

pdir 

K 7t/2 



if, 



if. 



(e")„ : V - V> 
(e% : i/j - i> 

{e w ) c : V - e 2ie iP (4.24) 
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We have also listed above the transformation for direct lattice rotations, which follows from 
the other results. Notice again that the transformation properties of if> under the space 
group are identical to the VBS observable. However, because of the non-zero charge of ip 
under U(l) c , we cannot generically identify if) with the VBS order. 
Finally, we note that the product of 7 mn and if) 

Pmn = lmn4>, (4.25) 

is indeed invariant under all the global U(l)'s, and so is the simplest composite operator 



which can serve as a density operator. From the space group transformations in Eqs. (4.22) 



and (4.24), we observed that p mn transforms as a linear combination of the components of 



the density at wavevectors (2irp/q)(m,n) + (ir, 0) and (2np/q)(m,n) + (0, 7r) 



The relationship in Eq. (4.25) is central to all our results: only the product of the 
monopole operator ip, and the vortex- anti- vortex composites represented by the 7 mn , is 
a physical observable. The requirement that we must consider the product of these dual 



operators can be traced to the constraint in Eq. (3.6) in the direct theory. There we noted 
that the monopole Berry phase was tied to the constraint on the sum of the spinon and 
doublon densities. In the dual theory if) = ipl4>l accounts from the spinon contribution, 
while 7 mn accounts for the density fluctuations in the paired doublon superfluid. 



2. Continuum theory 



We are now faced with the relatively straightforward task of writing down the most general 
action for the Cooper pair/doublon vortices, ip a e and the spinon vortices ipi,i, consistent with 



all the symmetries enumerated in Section IV A 1 



The quadratic kinetic terms are a direct transcription of the terms in Eq. (4.10), and lead 
to the Lagrangian 



Lq = \(dfj, - 2nia n - 2mc^)if)^\ 2 + {(d^ + 2ma u - 2mc lx )ipi\ 2 

+ |(<9 M - litibp + 2iric fJ )tp + t\ 2 + - 27ri6 M - 2iric fl )if-i\ 2 + .. . 



(4.26) 



Most crucial for our purposes will be the terms which directly couple the spinon vortices 
with the (fiai vortices. These are most directly deduced from Eqs. (4.22) and (4.24). Clearly, 
we need a combination of the 7 mn which transforms like the VBS operator under the space 
group operations so that the product with if) will be invariant under U(l) c and also under 
the space group. Such terms can only be constructed for even q, and were considered in 
Ref. EU in our present notation the simplest term is 



L 1 = X 1 



1p (TttO - ilOn) + H.C. 



(4.27) 
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Here we have labelled 7 mn by a subscript which identifies the associated wavevector 
(2Trp/q)(m,n), and will frequently use this notation below. A higher order term which 
will be important later is 



A 2 



^ (llo ~ ll) + H.c. 



(4.28) 



3. Phases 



We are now ready to use the vortex degrees of freedom to identify and characterize the 



phases introduced at the beginning of Section IV 



(i) AFM boson superfluid: Both the z a and the g± are condensed, and so all the vortex 
fields are gapped: 



(</> T > = , (^) = , =0 , = 0. 



(4.29) 



We will also need to consider independent condensates of bilinears of the vortices in the 
charges g± below, and so let us also note that in this phase 



(4.30) 



The low energy excitations of this phase consist of the 3 U(l) photons, a, b, c. These 3 
photons correspond to the 3 spin-wave modes that are easily deduced to be present in this 
phase of the direct theory. 

(ii) Paired boson superfluid: The restoration of the spin rotation invariance implies that the 
vortices in the spinons z a have condensed: 



(^t)^o , (Vi)^O , (<p +t )=0 , (<p-e) = 0. 



(4.31) 



The condensation of ^tj. implies that we also have (ip) ^ 0. However, this does not imply the 
appearance of VBS order, or broken translational symmetry, because of the non-zero U(l) c 



charge carried by ip. Note also that because of the coupling in Eq. (4.27), that a particular 
bilinear of the <p±e vortices has a non-zero condensate 



((ln0 ~ ilOir)) 7^ 0. 



(4.32) 



Again, this condensate does not break translational symmetry because it has to be combined 
with if} to obtain an observable neutral under U(l) c , and the combination is translationally 



invariant; indeed this translational invariance was used to derive the term in Eq. (4.27). All 



other linear combinations of bilinears of the ip±e vortices of the form in Eq. (4.30) have a 
vanishing expectation value. 

We can also use this vortex formulation to analyze the transition between the AFM boson 



29 



superfluid and the paired boson superfluid. Because the vortices <p±e are gapped in both 
phases, we can set tp±i = in all terms in the action. The critical theory then consists of 
2 complex scalars ^t.j. coupled to 2 U(l) gauge fields aj = a + c and a± = —a + c. The b 
gauge field is not Higgsed in either phase, and so remains gapless across the transition; this 
is just the Goldstone mode of the superfluid order, which is not connected with the critical 
theory. We can 'undualize' each complex scalar + U(l) gauge field combination: by the 
Dasgupta-Halperin dualitj^Sl this yields a critical theory of 2 complex scalars (and no gauge 
fields) with 0(2) x 0(2) symmetry. This critical theory is simply the easy-plane limit of the 
0(4) theory discussed in the direct formulation above and in Fig. |4j it is obtained from the 
models there by adding the easy-plane anisotropy term |2;-f| 2 |^| 2 . 

(Hi) AFM paired boson superfluid: Like the AFM boson superfluid above, all vortices have 
a vanishing condensate: 



<Vt) = o , (Vi) = o , {<p +t ) = Q , (M = o- 



(4.33) 



However, unlike the AFM boson superfluid, we should only allow for a condensate of the 
product g+g~, and not for the individual boson factors. In the dual variables, this means 
at least some linear combinations of the vortex bilinears (p_ m <p + £ should have a non-zero 
condensate while all the bilinears (p_ m (p + e have a vanishing expectation value. The simplest 
choice is to allow 

(Too) ^ (4.34) 

This does not break translational symmetry, and serves the important purpose of Higgsing 
out the c gauge field and ensuring that there is no single boson condensate. However 
it is also possible to choose other 7 mn to have non-zero expectation values. This option 



will be discussed further in Section IV B where we will show that in general such choices 
do break translational symmetry, and so lead to supersolid order (along with the AFM 
order already present here). However, there will be a number of other choices, distinct 



from Eq. (4.34), which do not break translational symmetry. All these choices lead to 



AFM boson superfluids which are identical in the sense of symmetry, but do have distinct 
'topological order' associated with the alignment of the gauge-dependent condensates of 
Lp J L m ip + £. We will provide a complete listing of all such inequivalent AFM paired boson 



superfluids in Section |IVB[ they are distinguished, in particular, by distinct universality 
classes of transitions involving the loss of AFM order. One particular choice that does not 
break any translational symmetry follows from our construction of the low energy theory 



in Section IV A 2 choose the condensate as in Eq. (4.32) and then all gauge invariant 



condensates that can be constructed out of it are translationally invariant. In this case, we 



observe from Eq. (4.27) that there is a term linear in the monopole operator ip in the action, 
and a consequent mixing of and 

(iv) Paired boson supersolid: Now we restore spin rotation invariance in the AFM paired 
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boson superfluid by condensing vortices in the spinons z a : 

(^T)^O , (^)^O , (<p+4)=Q , (<p- t ) = 0. (4.35) 
We assume the Cooper pair vortex bilinear condensates that are present here are the same as 



those in the AFM paired boson superfluid. For the choice as in Eq. (4.34), the present phase 
will have VBS order, as in the insulator (see Fig. [7]). This is because now the combination 
Poo — ipi^iloo has a non-zero expectation value, and this is a gauge-invariant observable 



which transforms like the VBS order parameter (see Eq. ( 4.25[ )). The transition to this VBS 



supersolid state from the AFM state in (Hi) will be of the easy-plane CP 1 variety, just as 
in the insulator. The second possibility noted above was to choose the Cooper pair vortex 



condensate as in Eq. (4.32): in this case the present phase will not break translational 
symmetry, and will indeed be identical to the paired boson superfluid in (ii). For the 
transition between the AFM paired boson superfluid and the paired boson superfluid (i.e. 
between the states in (iii) and (ii)), the mixing between ip^ and ?/>j implies that it is described 
by a single complex scalar (which is a linear combination of the these fields) coupled to the 
a gauge field. By Dasgupta-Halperin duality 58 , this transition is in the 0(2) universality 
class (and in the 0(3) class with full spin-rotation symmetry) — this is then nothing but 
the conventional SDW transition. We will consider other choices for the condensates of 



7 mn in Section IV B below, and also find one in which the paired boson superfluid is a 
nematic, for a case in which the density in the corresponding state in (Hi) has full square 
lattice symmetry. Other cases lead to a variety of supersolids with density oscillation periods 
which are characteristic of the total boson density 1 + x. 



B. Loss of AFM order in the paired boson superfluid 

This subsection deals with the nature of the superfluid phases (Hi) and (iv) above, and of 
the quantum phase transition between them. Both phases have a paired boson condensate, 
but no single boson condensate, and so the flux quantum is h/(2e). The first phase also has 
antiferromagnetic (AFM) order, and spin rotation invariance is restored in the transition 
to the second phase. We will see here that a rich variety of cases are possible, and we will 
present a few illustrative examples. 

We will argue that generically both the magnetically ordered and disordered phases break 
lattice symmetries. It is possible that the pattern of lattice symmetry breaking on the two 
sides of the magnetic phase transition is the same, in which case we expect a critical point 
in the 0(3) universality class; this includes the case where there is no lattice symmetry 
breaking in either state, as just noted above. It is also possible to have a phase transition 
where antiferromagnetism is lost, but a larger subgroup of the lattice symmetry is broken. 
Such transitions will either be first order or exotic (for instance, of a deconfined variety). 
We construct some specific examples of various scenarios. 
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We would like to describe the paired boson phase, in which (g+g~) 7^ 0. This operator 
corresponds to a monopole of the b field, hence b must not be Higgsed in this phase. On the 
other hand, we would like to suppress the single boson condensates z s g\ and z s g_, which 
correspond to monopoles of flux 1/2 in all three of the gauge fields. This can be achieved 
by Higgsing the gauge field c, which leaves the freedom for the gauge field a to be either 
in the Higgs or Coulomb phase - i.e. we can consider the loss of antiferromagnetism in the 
presence of a paired boson condensate. 

We will Higgs c by a condensate (e^ /3 + _/3 ~- 1 ) 7^ 0. Note that this condensate is not charged 
under the b field, and hence b remains unhiggsed as desired. If we, instead, condensed e*^+ 
and e ll3 ~ , independently, then both b and c would be Higgsed and the resulting state would 
be an insulator. 

We would now like to integrate the fluctuations of <p and c fields out, obtaining an effective 
theory for the spinon vortices = e*" 1 * and ipi = e lcti interacting with the gauge field a. 
In principle, the massless gauge field b, corresponding to the superfluid goldstone, also has 
to be included in the effective theory, however, as argued previously, it will decouple at low 
energies. The resulting theory will have two phases. In one phase, {ip-\) = 0, (ipi) = and 
a is massless - this is the antiferromagnetic phase. In the other phase, the spinon vortices 
condense, (ip-\) 7^ 0, (^4,) 7^ 0, the gauge field a is Higgsed and antiferromagnetism is lost. 

It is clear that for a generic set of condensates (7 mn ) the lattice symmetry is broken 
because we can construct gauge invariant observables like 7 mn T m i n i which transform non- 
trivially under the space group symmetry. This has important consequences for the structure 
of the effective theory for ip^, fields governing the loss of antiferromagnetism. The con- 
tinuum action for these fields will have the form, 



(4.36) 



L ° = 2P {V X + l( ^ " 2iria M 2 + \( 9 » + 2 ™ a M* + U(\^\ 2 , iVil 2 ) (4-37) 



where U is a potential term invariant under if)^ <-> V|- -^0 contains the lowest dimension 
operators invariant under independent phase rotations of ip^ and ipi^ However, due to the 



presence of the last term in (4.10), or more physically, due to the compactness of the direct 



gauge field A, only the combination (4.17) is a symmetry of the theory. Thus, we have an 
additional term L m , generated by the monopoles of the direct theory, which will break the 
"flux symmetry" 

17(1)* : Vt - e i0 ' 2 4> h ^ ^ e^Vi (4-38) 

The simplest terms in L m will be polynomials in the monopole operator ip = t/'t^i? which 
transforms as ip — > e ie ip under £7(1)$. In the well understood case of a pure spin system, 



lattice symmetry (4.24) implies that only "quadrupled monopoles" survive and the lowest 
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order term allowed in L m is, 

L m = ~y m {^ + (^) 4 ) (4.39) 



Intuitively, this quadrupling is due to the presence of oscillating Berry phases Q in (4.10) 



which lead to a destructive interference of single monopoles. We see that the action (4.39) 



preserves a Z 4 subgroup of £7(1)$, which from the transformation in Eq. (4.24) is iden- 



tified with (direct) lattice rotations. More generally, as noted in Section IV A, ip has the 
transformation properties of a valence-bond-solid (VBS) order parameter. On the magneti- 
cally ordered side of the phase diagram, the Z4 symmetry is unbroken, (-0) = 0, and hence 
the lattice symmetries are preserved, while the magnetically disordered side breaks the Z 4 



symmetry via, (ip) = (ip^ipi) 7^ leading to a VBS order. The term (4.39) is expected to 
be irrelevant at the critical point where the magnetic order is lost and hence monopoles are 
suppressed at the phase transition. Undualizing back to the direct theory, we obtain a model 
where the gauge-field A is non-compact and the resulting critical point is of a "deconfined" 
variety. In particular, a direct second order transition between the two phases is allowed. 

The above picture is still expected to hold in the present model at zero doping in the 
presence of a paired doublon superfluid, as the condensate (e^ /3+_/3 ~^) does not break any 
lattice symmetries. However, once we go to doublon superfluid states at finite doping and 
develop condensates (7 mn ) - we generically break lattice symmetry. Hence, lattice symmetry 
will be broken both in the magnetically ordered and disordered phases. Moreover, the 



monopole term in the spinon vortex action is no-longer constrained by Eq. (4.24), and single 



monopole terms will be generated from the coupling in Eq. (4.27), 

L m = -(y m ^ + y* m ^). (4.40) 
The fact that the monopoles are no longer quadrupled is roughly due to the spatially oscil- 



lating nature of the condensate (e^ /3+_/3 ~- ) ), which cancels the Berry phases in eq. (4.10). As 
we know, the direct gauge theory with monopoles allowed is equivalent to the 0(3) a- model 
(or its easy plane counterpart in the present case). Thus, we will have a phase transition 
in the 0(3) universality class (or 0(2) class in the easy plane case). This is consistent with 
our expectations since only the Neel order is lost and no new lattice order is gained across 
the phase transition. 

We note that though the above scenario is the most general one at finite doping, it is 
possible that the set of non-zero condensates (7 mn ) does not break the lattice symmetry, or 
break it only partially. Note that for a lattice operator g to be unbroken, it is enough that 
the product of g and a gauge rotation be preserved. So if the lattice symmetry is preserved 
by ('Jmn) U P to a rotation in the gauge group U(l) c , the transformations of the monopole 
field ip under the unbroken symmetries might be modified from the ones at zero doping in 



eq. (4.24). This unleashes a whole set of different possibilities for phase transitions out 



of the antiferromagnetic phase, accompanied by breaking of additional parts of the lattice 
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symmetry. We will present some examples of this scenario below. 

However, first we would like to discuss an alternative way of looking at the magnetically 
restored phase, where the spinon vortices ip^, ip^ are condensed. So far, we have been think- 
ing about the way the condensation of vortex-antivortex pairs (^_ i yj +n ) affects the spinon 
vortices. This is the correct logic for studying the phase transition where antiferromag- 
netism is lost. However, once the spinon vortices (or more generally the monopole field ip) 
are condensed, it is instructive to ask the reverse question: how are the Cooper pair vortices 
affected? The condensates if)^, ^ Higgs the c (and a) gauge fields and appear to break 
lattice symmetries. However, as long as K^t)! = K^J.)! (i- e - time reflection symmetry is 
unbroken), a combination of lattice and U(l) c rotations is always preserved. In particular, 
by a gauge rotation we can choose ip^ = ip± real, and then T x (i) c , T y , R ( ^f(e tn / 4: ) c , I x , T 
are preserved. Under these symmetries, the Cooper pair vortices transform as (we list only 
translations here for brevity), 

T x = T x (i) c : tp + i -> i<f- t i + i, <p-i -> -i(p +> i+i 

f y = T y : <f +i -> u~ l ip-i, <f-i -> Ld~ l ip + i (4.41) 

We observe, 

Tyf X = {U) b f x fy (4.42) 

q = - u = e 27rip/ « (4.43) 

V 1 + X /a a a\ 

- = (4.44) 

q 2 1 } 

Hence, the condensed monopoles endow the Cooper pair vortices with a new projective 



with, 



where 



implementation of the lattice symmetry. Comparing with eq. (4.13) we conclude that the 
condensed monopoles shift the effective density of boson from x to 1 + x, as we claimed 
above. 

Now, we come back to specific examples of phase transitions out of the antiferromagnet. 



In Sections IV B 1 IV B 2 and IV B 3 we will provide a complete listing of AFM paired boson 
superfluids which have full square lattice symmetry in the density; these cases will lead to 
transitions to the paired boson superfluid as shown in Fig. [5| The remaining Sections IV B 4 



IV B 5 will consider cases in which square lattice symmetry is broken in both the AFM and 
non-magnetic phases, but with the nature of the lattice symmetry breaking changing across 
the transitions. 
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(a) (b) 

FIG. 7: Orderings in the non-magnetic state in the presence of the (700) condensate (all q), or the 
(Ttttt) condensate (q/2 even). 

!• (Too) 7^ 0: Deconfined critical point to a valence-bond-solid (x - arbitrary) 

We imagine that the vortex-antivortex condensate present is (700)- As already noted in 



Section [IV A[ this condensate does not break any lattice symmetries. Hence, the discussion 
given for zero doping applies here. Namely, the antiferromagnetic state will break no lattice 
symmetries, and we will have a transition to a valence-bond-solid state via a deconfined 
quantum critical point described by the CP 1 field theory. The possible patterns of square 
lattice symmetry breaking are as in the insulator, and illustrated in Fig. [7} 



%■ (7vr0>, (70tt) 7^ (q - even) 

The simple example below illustrates the range of possibilities for the phase transition 
out of a paired doublon antiferromagnet. 

We assume that q is even and (7^0), (7ott) are non-zero. However, a particular choice 
°f (Ttto)) (7ott) exists where by combining lattice operations with U(l) c transformations, 
all lattice symmetries can be preserved. The resulting symmetries are T x (i) c , T y , i£ , T 
and i?^2 1 ( e± * 7r/ ' 4 )c ( we will drop transformations under the time reflection symmetry below, 
since it is never broken by (7 mn )). As we will see the factors (e ±Mr//4 ) c for corresponding 
to the presence of condensates ((7^0 =F ijon)) lead to two inequivalent scenarios. The new 
transformations for the monopole operator if), therefore, are 

■0 — > ^ 
1p —> 1p* 
>l{j -> ±1p ] 

i) -> v ( 4 - 45 ) 

Note the two different possible transformations under R^jf- If we choose the + sign, then 



T 
1 y 

odual 
rdual 
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(a) 



(b) 



FIG. 8: Ordering in the non-magnetic state in the presence of (7^0), (tott) condensates. 



the Berry phases of spinon vortices and Cooper pair vortices cancel each other, and the 
lowest allowed term in L m is, 



-y m (ip + ip 



(4.46) 



We can also view this term as arising from Eq. ( |4.27[ ). Thus, we obtain a theory with 
unsuppressed monopoles, and expect a phase transition in the 0(3) (0(2) in the present 
easy plane model) universality class. Note that the lattice symmetry will be unbroken on 
both sides of the phase transition. This is the conventional SDW transition between AFM 



paired boson superfluid and paired boson superfluid discussed in Section IV A 

Alternatively, if we choose the — sign for R^jf, the Berry phases of ^'s and ip's add-up 



and Eq. (4.45) become the transformations of a monopole operator in an antiferromagnet 
with odd-integer spins. The lowest allowed term in L m is, 



-y m {i? + (^) 2 ) 



(4.47) 



and the monopoles are "doubled"; this term arises from Eq. (4.28). The residual Z2 flux 



symmetry corresponds to direct lattice rotations. Thus, in the antiferromagnetic phase (ip) = 
and the lattice symmetry is unbroken. In the non-magnetic phase, (ip) 7^ and lattice 
symmetry is broken - there are two different patterns for this depending on microscopic 
details (similar to dimer and plaquette states of a VBS). In one case, (-0) = (V^) and the 
only broken symmetry is R^ 2 (broken to (-R^/ 2 ) 2 )- This is the nematic superconductor, and 
a schematic picture of this state is given in Fig. [8k. In the other case, (ip) = —{ift) and the 
lattice symmetry is broken to T x T y , T X T~ X , R^f, -T^ ual - A schematic picture of this state is 
given in Fig. [8)3. As for the nature of the phase transition in this case, it is expected that 
the doubled monopole is a relevant operator, which can lead to a direct first order phase 
transition. 
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3. (Ttttt) 7^ (q - even) 

This is another case in which the antiferromagnetic superconductor has no density wave 
order. The nature of the non-magnetic supersolid requires a separate analysis for q/2 even 
and odd. 



For q/2 even, we observe from Eqs. (4.21) and (4.22), in a gauge where (7^) is purely 



imaginary, the square lattice symmetry operations are T x , T y , I x ual , and R^f(i) c . The new 
transformations for the monopole operator if), therefore, are 



T 

rjdual 
K n/2 

rdual 



4 



—iip^ 



(4.48) 



It is now easily seen that this case is the same as the supersolid for (700) 7^ above, and 
the VBS order is as in Fig. [7j The transition between the AFM and non-magnetic states is 
described by the deconfined CP 1 theory. 



For q/2 odd, we find from Eqs. (4.21) and (4.22), again in a gauge where (7^) is 

purely imaginary, the square lattice symmetry operations of the antiferromagnet are T x , 

T< rdual 1 
7/ 1 J- rr 



■yi ± X 



) c , and R^f, an d therefore 



T 
T 

ndual 
H n/2 

rdual 



iijy 



(4.49) 



Again, only quadrupled monopoles as in Eq. (4.39 ) are allowed. Now in a non-magnetic state 
with (-0) 7^ 0, as above, the symmetry of the state is different for axg((ip)) = 0, tt/2, tt, 37r/2 
and arg((0}) = 7r/4, 37r/4, 5tc/4, 7ir/4. For the first case, the state with arg((0}) = 
preserves Tj, T y , (R^jf) 2 , and I X 1T , while for the second case, the state with arg((^)) = 7r/4 
preserves Tj, Ty, F&f, and I x ual T x T y . Unlike other states we have considered here, these 
states cannot be constructed purely out of modulations of the bond energy variable Qij = 
Si ■ Sj, where Si is the spin operator on site i of the direct square lattice. Instead, we also 
need a directed bond variable Py = (Si ■ Sj)(Sf — S]) which is a spin singlet observable 
obeying = —Pji- Note that because we are considering a doped system, the on-site spin 
fluctuates between Sf = 0,3/4, and so Py is not identically zero. The spatial modulations 
in these variables in the states for q/2 odd are shown in Fig. [9j As was the case for q/2 even, 
the transition between the AFM and non-magnetic states is described by the deconfined 
CP 1 theory. 
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FIG. 9: Ordering in the non-magnetic state in the presence of (7^) condensate for q/2 odd. The 

and (b) arg( (?/>)) = 7r/4. Both states are 4-fold 

Si ■ Sj, while the arrows 



states are defined by Eq. (4.49) for (a) arg((-0)) 
degenerate. The line patterns indicate the link energy variable Qij 



indicate the directed link variable Pij = (Si - Sj)(Sf 



Note that the arrows do not imply spin 



or charge currents (time reversal symmetry is preserved), and the state can be fully characterized 
by modulations in the charge density along the links. 



4- <7±*/2,o> ^0 (q = (mod 4)J 

Here we consider a transition in the background of condensates (7+^/2,0) f° r Q divisible by 
4. As we will see, one of the possible magnetically disordered states in this case is a charge 
density wave with period 4. Such a state is actually observed on the hole-doped side of the 
cuprate phase diagram close to the doping x — 1/8 (p — 1, q — 16p^. We note, however, 
that we reach such a state only from an AFM superconductor which already has density 
wave order as in Fig. [7^,. 

It turns out that once the condensates (7^/2,0) are present, T x is automatically bro- 
ken, as can be seen from the transformation properties of the gauge invariant observable 
(7 7r / 2 qi—k/2 0) ■ H° wever ; one can still arrange for T y and T%(i) c to be preserved (generally, 
the condensates ( r y±n/2,w) will also be allowed by these symmetries). Moreover, we have the 
choice of preserving rotations by 180° about either direct or dual lattice site. For brevity 
we will only discuss the later CctS6, clS it might be physically relevant. Then one can pre- 
serve T%(i) c , T y , (R^f) 2 and Z£ uaI . We recognize that a state with such a symmetry is a 
valence-bond-solid, see Fig. [7^, with superposed antiferromagnetic order. The transforma- 
tion properties of the monopole operator under the remaining symmetry group are, 



rp2 . 

x • 






T ■ 

y ■ 




-> ^ 


1 pduaK2 . 
\ K ir/2> ■ 






rdual . 

X 







(4.50) 



Thus, we again have the case of "doubled monopoles," eq. (4.47). Note, however, that 
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(a) (b) 

FIG. 10: An example of a phase transition out of an antiferromagnet with dimer order as in Fig. [7^, 
to the non-magnetic states (a) and (b) above. 

now the residual Z2 flux symmetry corresponds to translations by two lattice sites along 
the x direction. So the antiferromagnet has (ip) = and carries a VBS order. Once 
antiferromagnetism is lost and (xjj) 7^ 0, we break an additional subgroup of the lattice 
symmetry. There are again two cases: 

(a) (if>) = )• Then the remaining symmetry group is T*, T y , (R^f) 2 and l£ . This 
state is a bond centered charge-density-wave with period four. A cartoon picture of this 



state is shown in Fig. 10 1. Precisely such a configuration is observed by STM experiments 
on hole doped cuprates near x = 1/8.^ 

(b) (i/j) = —{ij}*). Then the remaining symmetry group is T 2 T y , TjT" 1 , (R^f) 2 and Z£ ual . 



A schematic picture of this state is shown in Fig. [TOp . 

Note that in both of the cases (a) and (b) we have a transition from a state with a unit 
cell area of 2 to a unit cell area of 4. 



5- (7±7r/2,o> + and (70^/2) / (q = (mod A)). 

The present example is of interest as we will be able to construct a deconfined critical 
point other than the one separating an antiferromagnet and a valence bond solid. 

Once (7±7r/2,o) an d (70^^/2) are turned on, both translations along x and y are broken. 
However, the combinations T 2 T 2 (i) c and T 2 T~ 2 (i) c are preserved. Moreover, one can now 
arrange for rotations about either direct or dual lattice site to be preserved. We will concen- 
trate on the direct lattice site case as it yields a deconfined critical point. Then, it turns out 
that we can mantain J^ ir and either (a) R^J 2 01 W •^^/2( e_w )c- The spatial modulations 



in such a state are shown schematically in Fig. 11 1. The transformations of the monopole 
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operator in the case (a) are, 



rplrpl . 

x y ' 






rplrp—2 . 

x y 
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+ -i> 


pdir . 
K ir/2 ■ 


^~ 




rdir . 

X 







with the case (6) differing only in the transformation under i?£y 2 > 



(4.51) 



(4.52) 



Therefore, "doubled" monopoles are permitted in case (b), making it of less theoretical 
interest. Below, we will, therefore, concentrate on the case (a). Here, as a consequence of 
a Z 4 flux symmetry associated with direct lattice rotations, only "quadrupled" monopoles 



are permitted as in Eq. (4.39), and a deconfmed phase transition may be possible. We note 
that we have addressed only the instability of the phase transition to monopole proliferation 
here. However, since the symmetry group in the present case is smaller than for the usual 
antiferromagnet - valence-bond-solid deconfmed critical point, other relevant operators can 
arise, such as e.g. (e T /^<9 M aj,)(|?/>f | 2 — IV'J.I 2 )- A classification and RG treatment of such 
operators is beyond the scope of this work. 

As for the pattern of spatial modulations: in the antiferromagnetic phase, we have (ip) 







and we obtain the state in Fig. 11 1. Once antiferromagnetism is lost and (ifj) 0, we break 
an additional lattice subgroup. As for the usual deconfmed critical point there are two 
possible cases: (ip) ~ 1 and (ip) ~ e Mr/ ' 4 . In the first case, the remaining symmetries are T, 



X 1 



T yi ( R f/2) 2T x T h J ? r ' alld We § et the State 111 Fi §- 
^4 
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d. In the second case, the remaining 



symmetries are T 4 , T 4 , (R^J 2 ) 2 T^Ty , I^ ir R^J 2 (for brevity we omit a figure of this state). 



V. DISCUSSION 



We have discussed different possibilities for the destruction of Neel order in metallic or 
superconducting two-dimensional quantum antiferromagnets by doping in a small density of 
charge carriers into the parent insulators. We have summarized our results already in detail 
in the introduction, and so we will be brief here. 

The standard SDW theory for the appearance of Neel order in a metal, generically requires 
an intermediate state between the large Fermi surface metal at overdoping and the small 
Fermi pocket state at very low doping. This intermediate state has 8 zero crossings in the 
fermion dispersion along the Brillouin zone diagonals. Because such an intermediate state 
has not so far been observed, we have examined other routes to connecting such states. We 
used a formalism which decomposed the electron operator as a product of bosonic spinons 
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(a) (b) 

FIG. 11: A candidate for a new deconfined phase transition from magnetic state (a) to a non- 
magnetic state (b) with a higher degree of lattice symmetry breaking. 

and fermionic spinless doublons. Despite our use of this 'fractionalized' approach, one of 
our results was the remarkable reappearance in this formalism of the conventional SDW 
criticality for the loss of Neel order in the superconductor. In addition, we also found other 
universality classes for the loss of magnetic order in the AFM+SC state which mimic those 
found in insulating antiferromagnets for different values of S, as shown in Fig. |5j For the 
metallic case, we found a transition to an algebraic charge liquid - the doublon metal. 

Our endeavour was motivated by the fairly strong evidence for a magnetic quantum 
critical point at which Neel order is lost in the electron-doped cuprates 1 1 * 2 * 3 * 4 * 5 * 6 !. We hope 
that the scenarios presented here will be tested in future experiments. A clear strategy to 



do so has been provided in Sec. I A 



In principle, our results here can be extended to the case of the hole-doped cuprates, 
which were considered earlier in Refs. IT%f4T)l The main phenomenological difficulty, as we 
noted in Section |TJ is that the antiferromagnetism in the La based hole-doped cuprates does 
not remain pinned at (tt, tt). However it may well be that (tt, tt) antiferromagnetism is more 
important in the other hole-doped cuprates. So, we can consider the transition from the 
antiferromagnetic metal to the holon metal — all of our analysis here on the transition to 
the doublon metal carries over, and the transition is in the 0(4) class. Unlike the doublon 
superconductor, the holon superconductor is not immediately unstable to confinement. The 
holon superconductor has Nf = 4 gapless Dirac fermion excitations which carry the U(l) 
gauge charge, and which suppress monopole proliferation for large Nf. It was assumed in 
Ref. HH1 that Nf = 4 was large enough for monopole suppression. However, in the event 
Nf > 4 fermions are required, the holon superconductor would be unstable to supersolid 
states, as discussed in the present paper. However, an understanding of the nature of the 
symmetry breaking in these supersolids requires computation of the monopole symmetry 
properties in the presence of gapless Dirac fermions — this we will address in future work. 
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APPENDIX A: TIME REVERSAL SYMMETRY 

Here we outline how the time reversal symmetry was implemented in the symmetry tables. 
We begin by defining time reversal on the lattice Grassman numbers c a , c°^\ 

T[c a ] = -eapeP (Al) 
T[c at ] = e a(3 cp (A2) 

This definition results in: 1) the dynamic term in the action is left invariant under time 
reversal, 2) the local electron density is invariant under time reversal, and 3) the electronic 
spin density changes sign under time reversal. 

Now using the following transformation for the bosons: T[b a ] = e a @b^ and T[6 Q ] = e a ^bo 
(for bosons the conjugates are of course not independent), we can infer what the g should 
transform into under time reversal [since we know how to write c in terms of g and b and 
we know how the c transform] . This is recorded in Table [IlJ The last step is to go from the 
lattice /, g and b into their continuum counter-parts and this requires knowledge of where the 



g fields have their minima, but is otherwise straight forward. This is recorded in Table III 
Note that unlike the analysis in Ref. HH1 the two sublattice fermions g± transform in the 
same way. This is related to the position of the doublon pockets in the BZ. 



APPENDIX B: MONOPOLE ACTION IN THE DOUBLON METAL 

This appendix will consider the action of a monopole in the U(l) gauge theory of the 



doublon metal state of Section III A The g± Fermi surfaces have low energy excitations which 
carry a U(l) gauge charge, and we will discuss their influence on the monopole dynamics. 
This problem was originally considered by Herbut et alW* using a duality analysis, but an 
oversight in their reasoning was pointed out by Hermele et alW^. Here, we will update the 
analysis of Herbut et al. and find that the action of a monopole diverges linearly with system 
size, consistent with other investigations^^. 
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We begin with the effective action of the U(l) gauge field, A^, after the g± fermions and 
the z a have been integrated out. At quadratic order, this can be written 56 in terms of the 
components of the 'electromagnetic' field F^ u = d^A u — d v A^\ 



Se S (A) 



d 2 kduj 



{e{k,u)Fl + yL{k,Lo)Fl y } 



(Bl) 



where e(k,u) and n(k,u>) are the dielectric constant and the magnetic permeability respec- 
tively. For \oj\ < Vpk, in the doublon metal we have®' 



(Jt(k,u) ~ 1 + x 



e(k,uj) 



1 (y ^ 'M 

¥ [ 1 + X T 



(B2) 



while at the 0(4) quantum critical point to the AFM metal the critical z a spinons lead to 



the propagator in Eq. (3.1) which corresponds to the magnetic permeability 

H(k,u) 



1 

k 



\U)\ 



X 



k 3 



(B3) 



We now apply the duality methods discussed in Section IV A to Eq. (Bl), and obtain the 
dual theory for the 'height' field h 



d 2 kdu 



k 2 



+ 



I? 



(B4) 



87r 3 \e(k,u>) fj,(k,u>) 

where the monopole operator is m ~ e 2mh . Note that for u — 0, the leading /c-dependence 
in Eq. (B4) is of order k 4 . Herbut et alW^ argued that renormalization effects would always 



generate an analytic term of order k 2 , and proceeded to investigate its consequences. As 
noted by Hermele et al^, this is incorrect — the leading term remains ~ k 4 because it is 
protected by the presence of the g± Fermi surface. 

We can now estimate the action of a monopole in the present Gaussian/RPA 
approximatio 



exp (— S m ) = (m) = exp (—2n 2 (h 2 )) 



from which we obtain 



2tt z 



d 2 kdu fi(k,uj)e(k,uj) 
8tt 3 k 2 fi(k,uj) + uj 2 e(k,Lu) 



(B5) 



(B6) 



From Eqs. ( |B2|B3[ ) we now find an infrared divergence ~ J d 2 kk 3 , which indicates that S m 
diverges linearly with system size. This justifies the neglect of monopoles in the doublon 
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metal and at the 0(4) quantum critical point. 
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